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To Vadim Schechtman, most cordially 

ABSTRACT. It is well known that the so called Bianchi IX spacetimes with 
SO(3)-symmetry in a neighbourhood of the Big Bang exhibit a chaotic behaviour of 
typical trajectories in the backward movement of time. This behaviour (Mixmaster 
Model of the Universe) can be encoded by the shift of two-sided continued fractions. 

Exactly the same shift encodes the sequences of intersections of hyperbolic geodesics 
with purely imaginary axis in the upper complex half-plane, that is geodesic flow on 
an appropriate modular surface. 

A physical interpretation of this coincidence was suggested in [MaMarlf]: namely, 
that Mixmaster chaos is an approximate description of the passage from a hot quan¬ 
tum Universe at the Big Bang moment to the cooling classical Universe. Here we 
discuss and elaborate this suggestion, looking at the Mixmaster Model from the per¬ 
spective of the second class of Bianchi IX spacetimes: those with SU(2)-symmetry 
(self-dual Einstein metrics). We also extend it to the more general context related 
to Painleve VI equations. 

1. Introduction, background and notation 

1.1. Plan of the paper. The Mixmaster Model of the early Universe with 
5'0(3)-symmetry in a neighbourhood of the Big Bang predicts a chaotic behaviour 
of “typical” trajectories (in the backward movement of time) encoded by the shift 
of two-sided continued fractions: cf. [KLKhShSiSS], [BoNo73], [May87], and refer¬ 
ences therein. 

The same shift encodes the sequence of intersections with purely imaginary axis 
of hyperbolic geodesics in the upper complex half-plane, see [Se85]. 
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This coincidence invites a closer attention, becanse the accepted mathematical 
sonrce of the classical Mixmaster chaos involves the behavionr of separatrices on 
the real bonndary of the respective dynamical system (cf. [Bo85]). Geometry of 
these separatrices and approximate dynamics that it encodes are not visibly related 
to hyperbolic geodesics. 

A physical interpretation of this coincidence was snggested in [MaMarld], Here 
we discnss and elaborate this snggestion, looking at the Mixmaster model from the 
perspective of Bianchi IX model with 5'17(2)-symmetry. 

More precisely, according to [MaMarld], the Mixmaster “classical chaos” shonld 
be considered as an approximation to an unknown quantum description of the tran¬ 
sition from the inhnitely hot qnantnm Universe at the moment of Big Bang to the 
cooling Universe gradnally htting a classical Einsteinian model. Time axis at the 
moment of Big Bang is pnrely imaginary, and it becomes real dnring the observable 
history of Universe. 

We argned that a mathematical model of snch a transition explaining Mixmaster 
chaos consists in inverse Wick rotation of time axes mediated by a move of time along 
random geodesics in the complex hyperbolic half-plane or rather its appropriate 
modnlar qnotient. This passage to the modnlar qnotient was critically important 
for onr argnment. It was snggested by two initially disjoint evidences. The hrst 
one was P. Tod’s remark that a conformal version of cosmological time in the 
Friedman-Robertson-Walker models has a natnral strnctnre of the elliptic integral 
(cf. [MaMarld], sec. 4.2). The second evidence was a well known formal coincidence 
of two encodings: of Kasner’s trajectories, on the one hand, and of hyperbolic 
geodesics with ideal ends, on the other hand. 

In this paper we develop and present fnrther details of this pictnre. Namely, we 
now look at snch a transition from the side of “gravitational instantons” that is, 
self-dnal Einstein spacetimes with 5'17(2)-symmetry. Many snch spacetimes have 
a natnral complexihcation, in particular, time axis can be extended to the complex 
half-plane, whereas the instantons themselves are dehned by restricting time to the 
imaginary semi-axis. 

Following the behaviour of the respective models along oriented geodesics in time 
connecting imaginary half-axis with real half-axis, we get the new aspect of the 
Mixmaster picture. This is the main content of this note. 

Structure of the paper. In the remaining part of section 1, we introduce some 
basic notation and constructions. 
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Sec. 2 compares (and shows a satisfying agreement) the seqnences of Kasner 
eras in the classical Mixmaster models with seqnences of geodesic distances between 
consecntive intersection of a geodesic with sides of the Farey tessellation. Finally, in 
sec. 3 and 4 we stndy an “instanton analogne” of the seqnences of Kasner solntions 
determining chaotic behavionr in the classical Mixmaster model. 

1.2. Continued fractions. We denote by Z, resp. Z+, the set of integers, 
resp. positive integers; Q, resp. R is the held of rational, resp. real nnmbers. For 
X G R, we pnt [x] := max {m G Z | m < x}. 

Irrational nnmbers x > 1 admit the canonical inhnite continned fraction repre¬ 
sentation ^ 

X = ;-T— =: [/co, /ci, ■ ■ ■ ], kg & Z+ (1-1) 

in which ko := [x], k\ = [l/x] etc. Notice that onr convention differs from that of 
[KLKhShSiSS]: their [/ci, /c 2 , ...] means onr [0, /ci, /c 2 , ■ ■ • 

1.3. Transformation T. The (partial) map T : [0,1]^ ^ [0,1]^ is dehned by 


T : {x,y) ^ 



’ y [l/x] 


( 1 . 2 ) 


If both coordinates {x^y) G [0,1]^ are irrational (the complement is a snbset of 
measnre zero), we have for nniqnely dehned kg G Z+: 


X = [0, /co, ki, /C2, ■ • y = [0, k-i, k-2 ,...]. 


Then 


1 r 1 


[0,/Ci,/C2,...], 


1 

y + [l/x] 


1 

ko + y 


[0,/co,/c-i,/c-2,---]- 


On this snbset, T is bijective and has invariant density 

dx dy 

ln2 ■ (1 -1- xy)‘^ 


(cf. [May87]). 



4 


Thus we may and will bijectively encode irrational pairs {x, y) G [0,1]^ by doubly 
infinite sequences 

{k) := [.. .k-2,k-i,ko,ki,k2,.. .],ki G Z+ 

in such a way that the map T above becomes the shift of such a sequence denoted 
T: 

T{k)s = ks+i- (1.3) 

1.4. Continued fractions and chaos in Einsteinian Bianchi IX type 
models. Bianchi classified metric space-times with a Lie group action transitive 
on space sections. In particular 4-dim Bianchi IX models of space-time can be 
of two types: with the symmetry group SO{3) or else SU{2). In the first case, 
metric has Minkowski’s signature, whereas in the second case it is Riemannian. In 
sec. 1 we survey the now classical results about chaotic behaviour in the 5'0(3)-case 
(Mixmaster Universe) and prepare ground for the treatment of 5'17(2)-models. Sec. 

2 and 3 are dedicated to the SU (2)-case. 

Consider the real circle defined in by equations 

Pa+Pb+Pc = l, pI+pI+pI^^- (1-4) 

Each point of this circle defines a 4-dimensional space-time with metric of Minkowski 
signature dt^ — a{t)dx‘^ — b(t)dy‘^ — c{t)dz^ with scaling factors a, b, c: 

a{t) = t^°‘, b{t) = c{t) — ,t > 0. 
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The main result of a series of physical papers dedicated to the Mixmaster Uni¬ 
verse can be roughly summarized as follows. 

A “typical” solution 7 of Einstein equations (vacuum, but also with various 
energy moment um tensors) with 5'0(3)-symmetry of the Bianchi IX type, followed 
from an arbitrary (small) value to > 0 in the reverse time direction t -t- 0 , 
oscillates close to a sequence of Kasner type solutions. (See subsection 2.3 below 
qualifying the use of adjective “typical” in this context). 

Somewhat more precisely, introduce the local logarithmic time O along this tra¬ 
jectory with inverted orientation. Its differential is dVl := -—, and the time itself 

abc 

is counted from an arbitrary but hxed moment. Then —)■ -|-cx) approximately as 

—logt as t —)■ -|- 0 , and we have the following picture. 

As = —log t —)■ -|-cx), a “typical” solution 7 of the Einstein equations determines 
a sequence of inhnitely increasing moments flo < fli 7 * ‘ ^ bln ^ • • • and a 

sequence of irrational real numbers Un G (1, +cxd), n = 0 ,1, 2 ,.... 

The time semi-interval [O^, bl„_|_i) is called the n-th Kasner era (for the trajec¬ 
tory 7 ). Within the n-th era, the evolution of a, b, c is approximately described by 
several consecutive Kasner’s formulas. Time intervals where scaling powers (pi) are 
(approximately) constant are called Kasner’s cycles. 

The evolution in the n-th era starts at time with a certain value u = Un> 1 
which determines the sequence of respective scaling powers during the hrst cycle 

(1.5): 

u 1 + u nil -|- u) 

Pi 1 ~ 2 ’ ~r~ — 2 ’ P3 ~r~ — 2 

1 + u + 1 + u + 1 + u + 


The next cycles inside the same era start with values n = Wn — 1, — 2,..., and 

scaling powers (1.5) corresponding to these numbers, rearranged corresponding to 
a bijection (1,2,3) —)■ (a, 6 , c) which is in turn identical to the previous one, or 
interchanges b and c (see [MaMar02]). 

After kn := [un] cycles inside the current era, a jump to the next era comes, with 
parameter 

1 

— i" • (1- 

Un Un\ 


Moreover, ensuing encoding of 7 ’s and respective sequences ("Uij’s by continued 
fractions ( 1 . 1 ) of real irrational numbers a: > 1 is bijective on the set of full measure. 
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Finally, when we want to inclnde into this pictnre also the seqnence of logarithmic 
times starting new eras, we natnrally pass to the two-sided continned fractions 
and the transformationn T. See some details in the next section. 

1.5. Doubly infinite sequences and modular geodesics. Let H := 

{z G C, Imz > 0} be the upper complex half-plane with its Poincare metric 
|(i 2 ;p/|Im 2 ;p. Denote also by H := H U{QU{cxd}} this half-plane completed with 
cusps. 

The vertical lines Re^ = n,n G Z, and semicircles in H connecting pairs of 
hnite cusps {j>//q') with pq' — p'q = ±1, cut H into the union of geodesic ideal 
triangles which is called the Farey tessellation. 

Following [Se85], consider the set of oriented geodesics /3’s in FI with ideal irra¬ 
tional endpoints in R. Let /3-oo, resp. /3oo be the initial, resp. the hnal point of /?. 
Let B be the set of such geodesics with /3-oo G (—1, 0), /3oo G (1, cxd). Put 

(d—oo [0, IvQj /C— 1 , k—2: • • - jj (doo [^1 5 ^2 •> ^3? ■■•jj ki ^ (1.7) 

and encode (5 by the doubly inhnite continued fraction 

[.. .k-2,k-i,ko,ki,k2,...]. ( 1 . 8 ) 

The geometric meaning of this encoding can be explained as follows. Consider the 
intersection point x = x{/3) of (3 with the imaginary semiaxis in H. Moving along 
(3 from X to /^oo, one will intersect an inhnite sequence of Farey triangles. Each 
triangle is entered through a side and left through another side, leaving the ideal 
intersection point (a cusp) of these sides either to the left, or to the right. Then 
the inhnite word in the alphabet {L, R} encoding the consecutive positions of these 
cusps wrt /? will be ... Similarly, moving from /3-oo to x, we will 

get the word (inhnite to the left) .. 

We can enrich the new notation ... R^°R^'^R^"^... (called cutting 

sequence of our geodesic in [Se85]) by inserting between the consecutive powers of 
L, R notations for the respective intersection points of (3 with the sides of Farey 
triangles. So xq := x = x{(3) will be put between R^° and and generally we 
can imagine the word 

... L^-^X-iR^°xqL^^xiR^^X2L^^x^R^^ ... ( 1 . 9 ) 
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We will essentially use this enrichment in the next section. 

2. Hyperbolic billiard, geodesic distance, and cosmological time 

2.1. Hyperbolic billiard. We will hrst present a version of the geometric 
description of geodesic flow: an equivalent dynamical system which is the triangular 
hyperbolic billiard with inhnitely distant corners (“pockets”). 

Here we use the term “hyperbolic” in order to indicate that sides (boards) of 
the billiard and trajectories of the ball (“particle”) are geodesics with respect to 
the hyperbolic metric of constant curvature —1 of the billiard table. This is not 
the standard meaning of the hyperbolicity in this context, where it usually refers 
to non-vanishing Lyapunov exponents. 

2.2. Proposition, a) All hyperbolic triangles of the Farcy tessellation of H are 
isomorphic as metric spaces. 

b) For any two closed triangles having a common side there exists unique metric 
isomorphism of them identical along this side. It inverts orientation induced by H. 
Starting with the basic triangle A with vertices {0, l,icxD} and consecutively using 
these identifications, one can unambiguously define the map b : iL —)■ A. 

c) Any oriented geodesic on H with irrational end-points in R is sent by the 
map b to a billiard ball trajectory on the table A never hitting corners. 

All this is essentially well known. 

It is also worth noticing that although all three sides of A are of inhnite length, 
this triangle is equilateral in the following sense: there exists a group of hyper¬ 
bolic isometries of A acting on vertices by arbitrary permutations. This group has 
a unique hxed point p := exp(7rf/3) in A, the centroid of A. 

In fact, this group is generated by two isometries: z ^ 1 — z~^ and symmetry 
with respect to the imaginary axis. 

Three hnite geodesics connecting the centre p with points i, 1+i, respectively, 
subdivide A into three geodesic quadrangles, each having one inhnite (cusp) corner. 
We will call these points centroids of the respective sides of A, and the geodesics 
(p, i) etc. medians of A. 

Each quadrangle is the fundamental domain for PSL{2, Z). 

2.3. Billiard encoding of oriented geodesics. Consider the hrst stretch of 
the geodesic (3 encoded by (1.9) that starts at the point xq in (0,icx)). If /cq = l,the 
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ball along /? reaches the opposite side (l,icxD) and gets reflected to the third side 
(0,1). If /co = 2, it reaches the opposite side, then retnrns to the initial side (0, ioo), 
and only afterwards gets reflected to ( 0 , 1 ). 

More generally, the ball always spends ko nnobstrncted stretches of its trajectory 
between ( 0 , icxD) and (l,icxD), bnt then is reflected to ( 0 , 1 ) either from (l,icxD) (if 
ko is odd), or from (0, ioo) (if ko is even). We can encode this seqnence of stretches 
by the formal word oo^° showing exactly how many times the ball is reflected “in 
the vicinity” of the pocket ioo, that is, does not cross any of the medians. 

A contemplation will convince the reader that this allows one to define an alter¬ 
native encoding of {3 by the donble infinite word in three letters , say a, b, c, serving 
as names of the vertices { 0 , 1 , icxD}. 

2.4. Kasner’s eras in logarithmic time and doubly infinite continued 
fractions. Now we will explain, how the donble infinite continned fractions enter 
the Mixmaster formalism when we want to mark the consecntive Kasner eras npon 
the f-axis, or rather npon the 0-axis, where O := —log J dt/abc 

In the process of constrnction, these continned fractions will also come with their 
enrichments, and the first new resnlt of this note will compare this enrichment with 
the one described by (1.9). 

We start with fixing a “typical” space-time 7 whose evolntion with t -t-0 nn- 
dergoes (approximately) a series of Kasner’s eras described by a continned fraction 
[ko,ki,k 2 , ■ ■ ■], where kg is the nnmber of Kasner’s cycles within s-th era [O^, Og+i). 
We have enriched this encoding by introdncing parameters Ug which determine the 
Kasner exponents within the first cycle of the era nnmber s by (1.5). A fnrther en¬ 
richment comes with pntting these eras on the 0-axis. According to [KLKhShSi85], 
[BoNo73], [Bo85], if one defines the seqnence of nnmbers 5g from the relations 

O 5 - 1-1 [1 Sgkg(^tig -(- l/'(n 5 })] 05 , 

then complete information abont these nnmbers can be encoded by the extension 
to the left of onr initial continned fraction: 


..., k-i, ko, ki, k2, ■ ■ ■] 


( 2 , 1 ) 


in snch a way that 


dg - xt/{xt + x^ ) 
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where 

xt ^ [0,ks,ks+i,...], xj = [0,ks-i,ks-2, ■ ■ ■]■ ( 2 . 2 ) 

2.5. Theorem. Let a “typical” Bianchi IX Mixmaster Universe be encoded by 
the double-sided sequence (2.1). Consider also the respective geodesic in H with its 
enriched encoding (1.9). 

Then we have “asymptotically” as s ^ oo, s G Z_|_.- 

S —1 

log Ll 2 s/Tlo ^ 2 ^ dist {X2r, X2r+l), (2.3) 

r=0 


where dist denotes the hyperbolic distance between the consecutive intersection points 
of the geodesic with sides of the Farey tesselation as in (1.9). 

Proof. According to the formulas (5.1) and (5.5) in [KLKhShSi85], and our 
notation (2.2), we have 


2s 2s 

log 02 ^i/^ 2 o ~ -^log (X+Xp) = ^ log ([/Cp_i,/Cp_2,/Cp-3,...])-[/jp,/jp+l,/jp+2, ■■•])• 

p=l 

(2.4) 

On the other hand, according to the formula (3.2.1) in [Se85], we have 
dist {xq,xi) = ^log{[ko,k-i,k- 2 , ■ ■ ■] ■ [/ci,/c2, ■ ■ • ] ■ [ki, ko, k-i,...] ■ [/c2,/ca, ■ • • ]) 
and hence, more generally, 

dist {X 2 r, X 2 r+l) = 

-log([/C2r, k2r-l: /^2r-2, . . .]-[/^2r+l, /^2r+2, . . .]-[/^2r+l, /^2r, /^2r-l, . . .]-[/^2r+2, /^ 2 r+ 3 , ■■•])• 

(2.5) 

Inserting (2.5) into the r.h.s. of (2.3), we will see that it agrees with the r.h.s. of 
(2.4). This completes the proof. 

The formula (2.3) justihes identihcation of distance measured along a geodesic 
with (doubly) logarithmic cosmological time in the next section. 
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During the stretch of time/geodesic length which such a geodesic spends in the 
vicinity of a vertex of A, the respective space-time in a certain sense can be ap¬ 
proximated by its degenerate version, corresponding to the vertex itself, and this 
justihes considering the respective segments of geodesics as the “instanton Kasner 
eras”. 


3. Mixmaster chaos in complex time and Painleve VI 

3.1. Painleve VI. Contrary to the separatrix approximation methods, the 
results about encoding of geodesics (3 with irrational ends and formulas for the dis¬ 
tances between consecutive cutting points are exact, but we did not yet introduce 
analogs of space-times hbered over geodesics as theirs time axes. We will do it 
in this section. The respective space-times are (complexihed) versions of Bianchi 
IX models with SU{2) (rather than 5'0(3)) action, the so called gravitational in- 
stantons. An important class of them is described by solutions of the Painleve VI 
equation corresponding to a particular point in the space of parameters of these 
equations: for us, the main references will be [To94], [Hi95], and [BaKo98]. 

However, the hyperbolic billiard’s picture of sec. 2 can be lifted to essentially ar¬ 
bitrary Painleve VI equations, and we will start this section with a brief explanation 
of the relevant formalism. 

Equations of the type Painleve VI form a four-parametric family. If the param¬ 
eters (a, /?, 7 , 5) are chosen, the corresponding equation for a function X (t) looks 
as follows: 


d'^x _ 1 /1 1 I \ uxV n 1 1 \dX 

~ 2\X^ X -1^ X -t) \ ^) ~[l^t-l^X-t)~^ 


X{X - l){x-t) 


In 1907, R. Fuchs has rewritten (3.1) in the form 


o t ^-1 

a + /3^ + 7(^_l)2+'5(YV7)5J' 


(3.1) 


i(l-f) 


X d‘^ ^ ^ d 1 


<X,Y) 


dx 


x{x — l){x — t) 


-ar + /3^+7(^_j)2+(.5 


1 

r (V -t)2 


(3.2) 
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Here he enhanced X := X(t) to (X,Y) := (X(t),Y(t)) treating the latter pair 
as a section P := {X(t),Y{t)) of the generic elliptic cnrve E = E{t) : Y^ = 
X{X — 1)(X — t). The section can be local and/or mnltivalned. 

In this form, the left hand side of (3.2) which we denote n{P) has a beantifnl 
property: it is a non-linear differential expression (additive differential character) 
in coordinates of P snch that fx{P + Q) = /u(P) +n{Q) where P + Q means addition 
of points of the generic elliptic cnrve E, with inhnite section as zero. In particnlar, 
//(Q) = 0 for points of hnite order. 

To see it, notice that the integral in the l.h.s. of (3.2) is additive modnlo periods 
of onr elliptic cnrve, considered as mnltivalned fnnctions of t. These periods are 
annihilated by the Ganss differential operator which is pnt before the integral sign 
in (3.2). 

The right hand side of (3.2) looks more mysterions. In order to clarify its mean¬ 
ing, notice that n{P) is dehned np to mnltiplication by an invertible fnnction of 
t. 

If we choose a differential of the hrst kind u on the generic curve and the symbol 
of the Picard-Fuchs operator of the second order annihilating periods of w, the 
character will be dehned uniquely. Moreover, it is functorial with respect to base 
changes (cf. [Ma96], sec. 0.2, 1.2, 1.3). In particular, if we pass to the analytic 
picture replacing the algebraic family of curves E(t) by the analytic one Ej- := 
C/(Z -|- Zr) T E H, and denote by ^ a hxed coordinate on C, then (3.1) and 
(3.2) can be equivalently written in the form 


where (cto,..., eta) := (a, -/?, 7, | - 5) and (Tq, Ti, T 2, T3) := (0, l,r,l + r), and 





(3.4) 


Moreover, we have 

=4(p(z.t) -ei(r))(p(z,T) -e2(T))(p(z,T) - e3(T)) 


(3.5) 
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where 

ei(T) = p(y,T), (3.6) 

so that Cl + 62 + 63 = 0 . 

The family Painleve VI was written in this form in [Ma06]. It was considerably 
generalised by K. Takasaki in [TaOl], in particnlar, he fonnd its versions for other 
families of Painleve eqnations. 

Now, any mnltivalned solntion z = z{t) of (3.3) dehnes a mnlti-section of the 
family which is a covering of H. In particnlar, if we can control its ramihcation 
and monodromy, then we may consider its behavior over geodesics with ideal ends 
in H and stndy the relevant statistical properties. The most accessible examples 
are algebraic solntions classihed in [BoaOS], [LiTyOS] and other works. 

However, here we will retnrn to Bianchi IX models, which according to [Hi98] 
correspond to the eqnation with parameters (a,/ 3 , 7 , h) = (|, —|, |), solvable in 
elliptic fnnctions. We will skip the beantifnl twistor geometry bridging Painleve VI 
and Bianchi IX and simply reprodnce the relevant resnlts from [To94] and [Hi95], 
somewhat reworked and simplihed in [BaKo98]. 

3.2. SU{2) Bianchi IX metric and scaling factors. Consider the SU{2) 
Bianchi IX model with metric of the form 


9 = F 


2 2 
7 2 ^2 

dfU, + 777o + TTTo T 




w? 



(3.7) 


Here fi is cosmological time, {aj) are 5'17(2)-invariant forms along space-sections 
with dai = aj A at for all cyclic permntations of (1,2,3), and F is a conformal 
factor. 

By analogy with the 5'0(3) case and metric dt‘^ — a(t)‘^dx‘^ — b{t)‘^dy‘^ — c{t)‘^dz‘^, 
we may and will treat Wi (as well as some natnral monomials in Wi and F) as 
SU ( 2 ) -scaling factors. 

However, contrary to the 5'0(3)-case, generic solntions of Einstein eqnations 
in the 5'17(2)-case can be written explicitly in terms of elliptic modnlar fnnctions, 
whereas their chaotic behavionr along geodesics in the complex half-plane of time is 
only a reflection of the chaotic behavionr of the respective billiard ball trajectories. 

3.3. Theta—functions with characteristics. Explicit formnlas in [BaKo98] 
nse the following basic f un ction of the complex argnments i/U G F, 2; G C, with 
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parameters (p, q) called theta-characteristics: 

q]{z,ifi) := E exp{—7r(m + + 2Tii{m + p){z + q)}. (3.8) 

m^Z 

It can be expressed through the theta-function with vanishing characteristics: 

■dip, q]{z, ip) = exp {—Trp'^p + 2mpq} ■ -d[0, 0](z + pip + q, ip). (3.9) 

All these functions satisfy classical automorphy identities with respect to the action 
of PGL{2,Z). 

3.4. Theorem. ([To94], [Hi95], [BaKo98].) Put 

il{p,q]:=‘ll{p,q](0,ifl) (3.10) 



(3,14) 
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or 


A> 0 , qGR, p E -+ iR. 

(B) Consider now a different system of scaling factors 

W[ := -h2^1ogi?2, := ^-K2^1ogi?3, 

p + qo dp p + Qo dp 

:= — + 2^1ogi?4, 

p -\- Qo dp 


(3.15) 


and 


F' ■- C(p + qof IVi'iyilVl 


(3.16) 


(3.17) 


where qq^C E Ti, C > 0. 

The metric (3.7) with these scaling factors for real p > 0 is real and satisfies the 
Einstein equations with vanishing cosmological constant. 

We will now consider valnes of G A C R in the vicinity of ioo bnt not 
necessarily lying on the imaginary axis. Since we are interested in the instanton 
analogs of Kasner’s solntions, we will collect basic facts abont asymptotics of scaling 
factors for ip -E- ioo. 

For brevity, we will call a nnmber r G R general, if r ^ Z U (1/2 + Z). 

For snch r, denote by (r) G (—1/2, 0)U(0,1/2) snch real nnmber that r+mo = (r) 
for a certain (nniqne) mo G Z. 

3.5. Theorem. The scaling factors of the Bianchi IX spaces listed in Theorem 
3.4 have the following asymptotics near p = +cx).- 

(i) For A = 0.- 

IT 1 

(3.18) 


W( ~ -^-, 

' 2 ’ 2 ^ ^ + 


(a) For A < 0 and general p: 

IFi ~ 7r(p) exp {7rf((p) — p)}, IF2 ~ ilFs, 

IF3 ~ -2'Ki {p + 1/2) ■ exp {tt* sgn {p)q} ■ exp{7rp(|(p)| - 1/2)}. 


(3.19) 
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(in) For A > 0, real q and p = 1/2 + ipo^po G R; 

-Wi Trpota.n{n{q - Pop)}W 2 ^-W 3 , 

W 3 ^ 2ttpo ■ {cosTT{q - pop))~^. (3.20) 


Comments. Theorem 3.5 shows that for general members of all solution families 
from [BaKo98], after eventual sign changes of some BA’s and outside of the pole 
singularities on the real time axis, we have asymptotically W 2 = IT3, Wi 7^ IT2. 

In the next section, we will show that such condition, when it is satished exactly 
rather than asymptotically, allows one to quantize the respective geometric picture 
in terms of Connes-Landi ([CoLaOl]. This gives additional substance to our vi¬ 
sion that chaotic Mixmaster evolution along hyperbolic geodesics reflects a certain 
“dequantization” of the hot quantum early Universe. 

Sign changes alluded to above are allowed, since Babich and Korotkin get their 
much simpler formulas by cleverly extracting square roots from expressions given 
in [Hi95]. Moreover, in the second version of their paper posted in arXiv in March 
2014, they corrected the signs of Wi and C (cf. their Lemma 2). For our purposes, 
this is not essential. 

Proof of Theorem 3.5. Directly from (3.9)~(3.11), we obtain: 

■^2 = ^ exp{-7r(m-S ~ 2 exp{-7r///4}, (3.21) 

m6Z 

■ds = exp{—~ 1 -|-2 exp{— tt//}, (3.22) 

m6Z 

114 = exp{—7rm^/u}(—1)™" ~ 1 — 2exp{— tt/u}. (3.23) 

m6Z 


Therefore 


TT 


4 ’ 


—log^ls ~ -27rexp{-7r;u}, 
dp 


—log'd4 ~ 27 rexp{—TT/uj. 
dp 


— \0gd2 
dp 


rsj 
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From this and (3.16), (3.17) one gets (3.18) for A = 0. 

Now consider the case A < 0, p general. 

Then from (3.8), (3.10), and (3.14) one gets 

q] = E exp{— 7r(^Tn + 27ri(^Tn H- p)q} ~ 

m6Z 


exp{27ri{p)q} ■ exp{-7T{pf p}, 


(3.24) 


becanse for general p, the leading term of q] corresponds to the nniqne valne of 
m for which (m + p)‘^ is minimal, that is, eqnals (p)^. 

Hence 

d 

—■dip, q] ~ 27rt {p)exp{2TTi{p)q} ■ exp{—7r{p)‘^p}. (3.25) 

oq 

Thns, from (3.12), and (3.21)-(3.25) we obtain 

■ 27rt (p)exp{27rt(p)((3' + 1/2)} ■ exp{—7r(p)^/u}x 


^d[p,q +1/2] i 




exp {—Tcip} ■ exp{—27ri{p)q} ■ exp{7r(p)^//} = —7r(p) exp {7ri{{p) — p)}. 
Fnrthermore, 

,,, _ i, , |;^[p+l/2,<,+ l/2] 

2 2 ^ * e'’'r"d\p,q] 

- -2 exp{—TT/u/d} ■27rt (p + 1/2) •exp{27rt(p + l/2)(q' + l/2)} - exp}—7r(p + l/2)^iu}x 

exp {—'Kip} ■ exp{—27rt(p)(3'} ■ exp{7r(p)^/u} ~ 

-27r {p + l/2)exp {m[{p + 1/2) - p- sgn{p)q]} ■ exp{7r//(|(p)| - 1/2)}. 


Notice that exponential terms were rewritten nsing the identity 


(p + 1/2) 


(P) 


1 

2 


sgn(p). 
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Similarly, 


T.. 1 , , ^^[P + 1/2,?] 

W3 := --^2^3— —^- 

2 ^[p, q] 


— ■ 2exp{—7r/u/4} ■ 27ri {p + 1/2) • exp{27ri(p + 1/2)5} ' exp{—7r(p + l/2)^iu}x 

exp{—27ri(p)5} • exp{7r(p)^/u} ~ 

-27ri(p + 1/2) ■ exp{7risgn(p)5} ■ exp{7r|U(|(p)| - 1/2)}. 


Comparing expressions for 11^2 and W 3 , one easily sees that W 2 = ihCs, where 
the exact sign can be expressed throngh p and q. 

For the conformal factor (3.13) we then get the following asymptotic: 

^ _ 2 W1W2W3 ^ 
ttA (^logr?[p,5])2 


(v + 1/2^2 

A (^)2 -exp {(p + 1/2) + (p) + 2 sgn (p) q} ■ exp{ 7 rp (2 |(p)| - 1)}. 

Finally, pass to the case A > 0. Pnt p = | + tpo, Po ^ R- We have again to 
locate hrst the leading terms as p ^ +cxd in 

'd[p,q] = E exp{— 7 r(m + p)^p + 27rt(m + p)^}, 

m^Z 

and also respective terms when p and/or q are shifted by 1 / 2 . Obvionsly, they 
correspond to the minimal valnes of Re (m + p)^, resp. Re (m + p + 1/2)2, 
m G Z. Since 

1 1 

Re (m + p)^ = (m + -)^ - Pq, Re (m + p + -)^ = (m + 1)^ - pp, 

in the hrst case there are two leading terms, for m = 0 and m = — 1 , and in the 
second case jnst one, for m = — 1 . 

Thns, for A > 0, we have 

-d[p, q] ~ exp{7rp(po - 1/4)} ■ [exp{27rtp5 - Tvipop} + exp{27rt(p - 1)5 + Tripop}]. 
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The sum of two terms in square brackets can be rewritten so that in the end we 
obtain 


~ exp{7in{pl - 1/4)} • exp{-27rpo?} ■ cos7r{q - pop). (3.26) 

'd[p + 1/2, q] ~ exp{nppl} ■ exp{-2npoq}. (3.27) 

When we have to replace a real g by g + 1/2, we may do it formally in the right 
hand side expressions in (3.26), (3.27). 

Therefore, we have from (3.12), (3.22) and (3.27): 

W2__. V2,Q + 1/2] 

W 3 ^ ^ e^^P^i9lp+l/2,q] 


Now, 


exp{-27rpo(Q + 1/2)} 
exp{7ri(l/2 + ipo)} ■ exp{-27rpoq} 


W: 


3 •- 


1 a a ^"^^ + 1/2^] 

-l/ol/s — --;- 

2 ^[p, q] 


rsj 


exp{- 7 r/// 4 } ■ (27rpo) ■ exp{7rppl} ■ exp{-27rpoQ} 
exp{7rp{pl - 1/4)} ■ exp {-27rpoq} ■ cos 7 r(Q - pop) 


27TP0 ■ {cos7T{q - Pop)) ^ 


Furthermore, 


-IFi 


^ ^^bw + 1 / 2 ] 

2^4 


1 exp{7ip{pl - 1/4)} ■ ^[exp{-27rpo(Q + 1/2)} ■ cos7r(Q + 1/2 -poh)] 

2 exp{7r(i/2 -po)} ■ exp{7rp(pQ - 1/4)} ■ exp {-27rpoq} ■ cos7r(ij - poh) 

1 ^[exp {-27rpo(? + 1/2)} ■ cosTi{q + 1/2 - pop)] 

— • - rs^ 

2 exp {7r(i/2 — po)} ■ exp {—27rpoq'} ■ cos 71(5 — pop) 

1 ^[exp{-27rpo(Q + l/2)} ■sin7r(Q-pop)] 

2 exp {-27rpo(q'+ 1/2)} ■ cos 77(5 - Pop) 


rsj 
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7rpo tan{7r («? - po/i)} - 2 • 

This completes the proof of Theorem 3.5. 


4. Theta deformations of gravitational instantons 


4.1. Theta deformations. In Section 5 of [MaMarl4] we showed that the 
glning of space-times across the singnlarity nsing an algebro-geometric blownp can 
be made compatible with the idea of spacetime coordinates becoming noncommn- 
tative in a neighborhood of the initial singnlarity where qnantnm gravity effects 
begin to dominate. 

This compatibility is described there in terms of Connes-Landi theta deforma¬ 
tions ([CoLaOl]) and Cirio-Landi-Szabo toric deformations ([CiLaSzl3]) of Grass- 
mannians. 


Here we consider the same problem in the case of the Bianchi IX models with 
SU ( 2 )-symmetry, namely whether they can be made compatible with the hypothesis 
of noncommntativity at the Planck scale, nsing isospectral theta deformations. 


The metrics on the sections, in this case, are only left 5'17(2)-invariant. We 
show that among all the SU{2) Bianchi IX spacetime, the only ones that admit 
isospectral theta-deformations of their spatial 5'^-sections are those where the met¬ 
ric tensor 


g = W 1 W 2 W 3 dg? + 


W2W3 

Wi 


(jf -F 


WiWs 2 , 

- (To + 

W 2 


W1W2 

W 3 


(Tq 


(4,1) 


is of the special form satisfying wi ^ W 2 = W 3 (the two directions (T 2 and <73 
have eqnal magnitnde). In these metrics, the sections are Berger spheres. This 
class inclndes the general Tanb-NUT family ([TanbSl], [NUT63]), and the Egnchi- 
Hanson metrics ([EgHa79a], [EgHa79b]). The theta-deformations are obtained, as 
in the case of the deformations Sg of [CoLaOl] of the ronnd 3-sphere, by deforming 
all the tori of the Hopf hbration to noncommntative tori. 


4.2. Proposition. A Bianchi IX Euclidean spacetime X with SU{2)-symmetry 
admits a noncommutative theta-deformation Xg, obtained by deforming the tori of 
the Hopf fibration of each spacial section to noncommutative tori, if and only if 
its metric has the SU{2) x U(1)-symmetric form 


g = W 1 W 3 dg^ 


wt 

+ — 
Wi 


+ Wi (a\ + ( t |) 


(4.2) 



20 


Proof. In appropriate local coordinates the 5'17(2)-invariant forms {ai) satis¬ 
fying relations dai = aj A ak for all cyclic permntations {i,j,k) have the explicit 
form 

cTi = xi dx2 — X2 dxi -t- 2:3 dxQ — xq dxs = -{d'il: -|- cos 9 dcf)), 


02 — X2 dxs — x^dx2 + 2:1 dxQ — 2:0 dxi = - (sin '0 (i 6 * — sin 6 * cos '0 
Os — xs dxi — 2:1 dxs + X2 dxQ — 2:0 dx2 = cos tf) d 9 — sin 6 * sin '0 


with Enler angles O<0<7r, 0<(/)<27r and 0 < '0 < dyr (for the SU(2) case). 
The Hopf coordinates (^ 1 ,^ 2 ,^) are dehned by 


zi 2:1 -t- ix 2 — cos - = cosry, 

Q 

Z 2 Xs + ixQ — sin - = sinry. 

2 

Eqnivalently, identifying with nnit qnaternions, we write q G SU (2) as 


q := 




/ cos ?7 sin ?7 

y sin rj cos rj 


where \zi\^ -|- lz 2 j^ — 1 and (^i,^2,h) are the Hopf coordinates as above. 

The noncommntative 0-deformations ([CoLaOl]) of the 3-sphere are obtained 
by deforming all the 2-tori of the Hopf hbration to noncommntative 2-tori Tg. 
Namely, replace q with 

U COST] V sin ?7 
—V* sin ?7 U* COST] 

where U, V are the generators of the noncommntative 2-torns Tg. 

Then one obtains the algebra generated hj a = U cos ry and (3 = V sin rj with 
aj3 = a*(3 — a*a = act*, (3*j3 — PP* and aa* -t- PP* = 1. 

It is shown in [CoLaOl] that this deformation is isospectral with respect to the bi¬ 
invariant ronnd metric on in the sense that the data of the Hilbert space of 
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square integrable spinors H = L‘^{S^, S) and the Dirac operator D for the round 
metric on give rise to spectral triples on the deformed algebras Sg. 

In fact, the general result of [CoLaOl] shows that isospectral theta-deformations 
can be constructed whenever there is an isometric torus action. In particular, in 
our case the question reduces to whether the action of that rotates the tori of 
the Hopf hbration preserves the Bianchi IX metric. 

In Hopf coordinates the action of is given by (^ 1 ,^ 2 ) : (^ 1 ,^ 2 ) t + 

^ 2 ), or in terms of the Euler angles, {u,v) : (</>, e-)- (cf) + u,^|^ + v), with ti = 
{u + v)/2 and t 2 — {v — u)/2. It is immediate to check that the 17(l)-action 
u : (j) ^ + u leaves the 1-forms invariant. This is the 17(l)-action of the Hopf 

hbration ^ ^ S'^. The form ai is also invariant under the other 17(l)-action 

V : 'ip ^ 'ip + V, while the other forms a2, (Ts transform as 

v*a 2 = - (sin('0 + /3) d9 — cos{'ip + /?) sin 9 d(p) 

2 

v*a 3 = -(— cos{^p + (3) d9 — sin('0 + (3) sind d(p)^ 
hence it is clear that we have v*g = g for a Bianchi IX metric 

g = dg^ + af + cr| -|- (4-3) 


if and only if 6 = c. In the case b = c, with 

Qp" (P‘ 

g = dg^ -I- — {dip -)- cos 9 d(p)‘^ -I- -^{d9‘^ -I- sin^ 9 dcp^), 

the action is isometric and the resulting theta-deformations are therefore isospec¬ 
tral, with spectral triples {A, H, D), with A = C'°°(5'|), and spinors H = L‘^{S^, S) 
and Dirac operator D with respect to the Bianchi IX metric with b = c. 

This is in stark contrast with the situation described in [EsMarl3], where (Lorentzian 
and Euclidean) Mixmaster cosmologies of the form 

=F(it^ -|- a{t)^dx^ + b{t)‘^dy‘^ -|- c{t)^dz^ 

were considered, with T^-spatial sections, which always admit isospectral theta- 
deformations. 
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We have recalled in the previous section how the self-duality equations for the 
SU (2) Bianchi IX models can be described in terms of Painleve VI equations [To94], 
[Hi95], [Ok98], and how the general solutions (with wi ^ W 2 ^ ws) can be written 
explicitly in terms of theta constants [BaKo98], and are obtained from a Darboux- 
Halphen type system [PeVal2], [Tak92]. In the case of the family of Bianchi IX 
models with SU{2) x 17(l)-symmetry, considered in Proposition 42, this system has 
algebraic solutions that give 


1 

W2=W2, = -, Wi 

n- Ho 


H- H* 


(4.4) 


with singularities at (curvature singularity), //q (Taubian inhnity) and oo (nut). 
The condition < hq avoids naked singularities, by hiding the curvature singular¬ 
ity at behind the Taubian inhnity, see the discussion in Section 5.2 of [PeVal2]. 

Consider the operator 


Db ~ —i 


jXi 

Xo - iX. 


X2 + iXo 


+ 


A2 +2 
2A ^ 


(4.5) 


where {Xi, X 2 , X^} constitute a basis of the Lie algebra orthogonal for the bi¬ 
invariant metric. Assume moreover that the left-invariant metric on is diagonal 
in this basis, with eigenvalues {w‘^/wi, wi, rci}, with w = W 2 = ws and A = w/wi, 
and where the Wi are as in (4.4). Consider also the operator 


D 


1/2 

w^ w 




d 1 ,w Iwi. 

-h -(-^- - 

dfi 2 w 2wi‘ 


+ Wi Db\\=-ML- 

' W 2 


(4.6) 


4.2. Proposition. The operators D of (4-6) give Dirac operators for isospec- 
tral theta deformations of the SU{2) x U{1)-symmetric spacetimes of Proposition 

4 . 2 . 

Proof. We consider the frame 9^ with i G {0,1, 2, 3}, given by 

9^ = uw dpi, 9^=u\ai, 9‘^=ua2, 9^=uao, 

1/2 

where u = w-^ and \ = w/wi, ior w = W 2 = w^- Since the Oi satisfy dai = aj A ak 
for cyclic permulations {i,j, k} of {1, 2, 3}, we have d9^ = 0, and furthermore 

= iuX + u\) dp ^ ai + u\a 2 ^ (To = —(-+ ^) A ^ X^9^ A 9^, 

uw u X uX 
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d9^ = iidfx A a2 + uas A ai = — —9^ A9^ -- 9^ A 9^, 

uw u uX 

\ u 1 

d9^ = udn A as + uai A a2 = - 9^ A9^ -\ - - 9^ A 9^ 

uw u uX 

where dot denotes the time derivative. 

Proceeding then as in [ChCol2], we nse the d9^ to write the spin connection and 
we obtain a Dirac operator of the form 


D = -f^ 


1/2 
W, W 


d 1 ,w 1 rfii . 

-h -(-^- - 

dfi 2 w 2 wi' 


+ 


w 


1/2 


D 


w 


' wi 


or eqnivalently of the form (4.6), where is the Dirac operator on a Berger 
3-sphere. The explicit form of Dirac operator on a Berger 3-sphere with metric 
A^cr^ + ct| + (jg was compnted in [Hi74], and it is given by the operator (4.5). 

As in [EsMarl3], the Dirac operator of Proposition 4.3 can be seen as involving 
an anisotropic Hnbble parameter H. In the case of the metrics (4.3) of [EsMar]13 
this was of the form 


H 


l a b c 

o “ + I + “ 

3 \ a b c 


with a, 6, c the scaling factors in (4.3). 

In the case of the SU{2) Bianchi IX models, the anisotropic Hnbble parameter 
is again of the form H = ^(i^i + H 2 + i^s), where now the Hi correspond to the 
three directions of the vectors dnal to the 5'17(2)-forms ai in (4.1). For a metric of 
the form (4.2), or eqnivalently 


i2| 2\22| 22, 22 

g = uw dg + u X a^ + u a 2 + u as, 

with u, A, w as in Proposition 4.2, we take the anisotropic Hnbble parameter to be 

3\ uX u 6 \ u X 


ii Iwi X w wi 
u 2 wi ’ A w wi’ 


where 
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SO that 


H 


1 / w 1 \ 

3 2 / ’ 


as in (4.6), so that we can write the 4-diniensional Dirac operator in the form 


D 


-i"-^a4- + Ih]+Db, 


UW 


dfi 


where Db = /w) Db\x=:^ is the Dirac operator on the spatial sections 

with the left 5'17(2)-invariant metric. 

Notice that in the constrnction above we have considered the same modnlns 9 
for the noncommntative deformation of all the spatial sections of the Bianchi IX 
spacetime, bnt one conld also consider a more general sitnation where the parameter 
9 of the deformation is itself a f un ction of the cosmological time fi. 

This wonld allow the dependence of the noncommntativity parameter 9 on the 
energy scale (or on the cosmological timeline), with ^ = 0 away from the singnlarity 
where classical gravity dominates and noncommntativity only appearing near the 
singnlarity. Since a non-constant, continnonsly varying parameter 9 crosses rational 
and irrational valnes, this wonld give rise to a Hofstadter bntterfly type pictnre, with 
both commntativity (np to Morita eqnivalence, as in the rational noncommntative 
tori) and trne noncommntativity (irrational noncommntative tori). 

Another interesting aspect of these noncommntative deformations is the fact 
that, when we consider a geodesic in the npper half plane encoding Kasner eras 
in a mixmaster dynamics, the points along the geodesic also determine a family of 
complex strnctnres on the noncommntative tori T| of the theta-deformation of the 
respective spatial section. 
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